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ABSTRACT 

Comparison of two-phase choked flows in normalized coordinates were 
made between pure components and available data using a reference fluid 
to compute the thermophysical properties. The results are favorable. So- 
lution of the governing equations for two LNG mixtures show some possible 
similarities between the normalized choked flows of the two mixtures, but 
the departures from the pure component locii are significant. 
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NOMENCLATURE 

* 2.5 constants, *q. (11) 
b 1>4 constants, *q. (12) 

G mass flux (G * pU) 

G* mass flux normalizing parameter (|^P c e c /Z c ) 

G* for corresponding states with refer- 

ence fluid (0) and working fluid (a) 

3 configurational free energy 

h scaling function - volume 

f scaling function - temperature 

m molecular weight 

P pressure 

R gas constant 

S entropy 

T temperature 

U fluid velocity 

V specific volvne 

x mole fraction 

Z compressibility 


Greek : 

p density 

o,d binary Interaction parameter 

»,♦ configuration function 
« PI tier's eccentric factor 

i property parameter 

Subscripts: 
c critical 

e exit 

a fluld-a 

1j fluid 1 - fluid i 

e reference fluid 

m mixture 

t throat 

R reduced 

Superscript: 
c critical 

INTRODUCTION 

Using conventional corresponding states tech- 
niques, we have found normalizing parameters which 
can be applied in two-phase choked flow computations 
(refs. 1 to 4). The pure component properties of 
the working fluid, along with the flow equations, 
are used to define the critical mass flux ana pres- 
sure ratio profiles. When compared to available 
experimental data the agreement was good (refs. 4 
and 5). 

Although the normalized charts are applicable 
to a wide variety of simple fluids, it has not been 
demonstrated, via direct calculation, that a refer- 
ence fluid (0) can be used to predict the proper- 
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ties of Another fluid (•) In conjunction with calcu- 
1 At Ions of two-phased choked flow parameters for 
fluid (•) using fluid (0); the pure-pure problem. 

further, the eppllcetlon of these pAremeters to 
fluid mixtures end the problems essoclAted with the 
computetlon of two-phase choked flows of mixtures 
hes not been treeteo using the one fluid ven der 
Weals ApproximAtlon (refs. 6 to 11) for preoicting 
thermophysical properties, the pure-mix problem. 

In this psper, we wish to Address both the 
pure-pure end the pure-mix problems In two-phese 
choked flows using the thermophysic el properties 
cooes of references 12 end 10, respectively. It Is 
importent to study the fluid mechenics of mixtures 
as most working fluids ere mixtures of verying com- 
position. 

GOVERNING EQUATIONS 

The equetions governing homogeneous equilibrium 
two-phese choked flows ere oeteiled in references 2 
end 4. The mess flux is given as equetion (1): 



Constraints: 


defined in terms of Assumed mixing rules end enti- 
ce 1 peremeters releting component-k to the reference 
fluid (0): 



with binery interaction peremeters njj end 
Configuration functions • end * ere slowly vary 
ing functions of Vr end Tr end can be defined in 
a variety of ways. The most successful ere the leach- 
Leyland forms (ref. 6), or their modifications such as 
given by Mollerup (ref. 10) end McCarty (ref. 11) 

• • 1 ♦ (*if - M Q)£*j(Vp “ ” *3^R “ *4^ n ^R 

(ID 



S(Pq,Tq) • constant (3) 

The solution of equation (1) is easy enough if 

V is known. The whole problem lies in computing 

V subject to constraints 2 an: 3. For the simpli- 
fied pure-pure problem, this is rather straight- 
forward using a thermophysical properties package 
such as GASP (ref. 12), but can become complex for 
the pure-mix problem, depending on mixture composi- 
tion. 

In this paper, a reference fluid (0) will be 
useo to determine the thermophysical properties 
necessary to satisfy equations (1), (2), and (3) for 
both the pure-pure and the pure-mix problems. The 
mixture equations, detailed by Mollerup (refs. 7 to 
10), are given here in abridged form where l is 
the compressiblity function, $ the configurational 
Gibbs free energy, and x, the mole fraction of 
component j. 

ZJV.T.Xj) - 2 0 (V/h,T/f) (4) 

9 m* V,T,x j) “ ^(V/h.T/f) 

♦ RT^irf Xj In x j — In ^ (5) 

where the scaling functions, f and h are defined 
as: 

h • ££ Vj h ij < 6 > 

fh . EE VjV ij w 

* J 

and the interaction parameters hjj and ffj are 



The coefficients a, and bj are estimated from 
experimental data (refs. 10 and 11). 

For the pure-pure problem, equations (4) to 
(10) are simplified since the crossterms do not en- 
ter; for example, with only one component, equations 
(5) to (7) become: 


fl(V,T,l) - ffl 0 (V/h,T/f) - RT In h 

(13) 

h • . 

(14) 

' • (“X)‘ 

(15) 

For the derivatives necessary to solve 
erning equations, see Mollerup (ref. 10). 

the gov- 


RESULTS 


The pure-pure problem results are presenteo 
first since they are direct and readily explainec in 
terms of data. By contrast, the pure-mix problem 
which is complex, is further complicated by lack of 

data. 

Pure -Pure 

Using the thermophysical property codes of ref- 
erence 12, the governing equations were solved for 
pure components nitrogen ana oxygen on three selec- 
ted Isotherms for which accurate experimental data 
are available. Subsequent verification checks were 
made using the Mollerup code, (ref. 10), with simi- 
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lar results; the exception being near the saturation 
locus where some discrepancies noted. 

In figures 1 (a) end (b) the calculated critical 
mass flux and pressure ratio profiles in normalized 
coordinates are compared to experimental data. The 
curve labeled CH 4 - N? was calculated using the 
properties of methane, that is, fluid ( 0 ) is methane 
and fluid (•) is nitrogen; similarly, for the curve 
marked tHa - 0?. The use of nitrogen or oxygen 
as the reference fluid oio not alter the results of 
figures 1 to 3. The agreement with data is quite 
gooa, and similar departures from a single curve 
were found when using fluid nitrogen or fluid oxygen 
properties directly. 

Approaching the critical Isotherm, while the 
mass flux agreement is good, figure 2 (a), a rather 
large amt as yet unexplained departure in the pres- 
sure profiles between data ano theory can be noted, 
figure 2 (b). 

Close to the critical isotherm, the agreement 
with data is again good, as noted in figures 3(a) 
and (b). 

it must be pointed out that the fluid nitrogen 
and fluid oxygen data are in very good agreement. 

The Discrepancies lie in the governing equations, 
and their solutions. 

Pure-Mix 

Some initial computations using the governing 
equations and the more conventional corresponcing 
states approach to defining the normalizing param- 
eters appeared successful, but a closer examination 
revealed several problems ana the approach was aban- 
doned. Calculations for two mixtures are presented 
as a method of comparison, 

Conventional corresponding states techniques, 
some of which are given in reference 13, cannot be 
applied to mixtures with any degree of confidence. 
The linear molar fraction method; 


£ C - EtjXj (161 

i ' 

P c » Z C RT C /V C (17) 

where t c represents a property, is simple to 
apply but not reliable and the complex techniques 
recommended by reference 14 are better. 

With the molar fractions of mixtures Mixj and 
Mix? given in tables 1 and 11, some estimates of 
the critical parameters can be made which then can 
be compared to those estimated using Mollerup's pro- 
gram (ref. 10). The results are given in tables 1 
and 11 to illustrate the point. Further, the cir- 
conaenbar (maximum pressure) and circondentherm 
(maximum temperature) are given for each mixture. 

For all calculations the reference fluid is methane. 

The nature of the vapor pressure locus at a 
fixed composition is given as figure 4 for mixtures 
M1xj and Mix?, respectively. While alike in 
general appearance, the critical point of Mix? is 
close to the circondenbar. Further note that at 
temperatures greater than critical, following an 
isotherm from the high-pressure gas side with de- 
creasing pressure, condensation first occurs on the 
dew line and attains a maximum density followed by 
an evaporation at the dew line again, ano returns to 
the low pressure gas state. The composition of the 
fluid changes, as illustrated in figure 5 for the 
major components of Mix}. At lower pressures, the 


liquid is dominated by the least volatile components 
n-heptane and hexane. Finally, a temperature en- 
tropy diagram, figure 6 , for M1xj can be useo to 
"spot" ano follow the expansion process. Note that 
expansions occurring between 190 ano 400 K, depend- 
ing on pressure, terminate at about the saturation 
locus near the critical point to the circonaentherm; 
that is, the region most difficult to handle and 
which is the one of most interest. 

With this behavior, it should not be too 
surprising that we have significant difficulties in 
computing the critical mass flux for Mix^ and Mix?. 

The critical mass flux for Mixi in normalized 
coordinates is presented in figure / with the con- 
ventional pure-pure normalized locii, as back- 
ground. There are rather significant differences: 
J \) the saturation locus, ( 2 ) the slopes of the iso- 
therms, and (3) regions where computations are in- 
complete; the exception being the gaseous phase 
where similar behavior is noted. In many computa- 
tional attempts, no reasonable solution was found as 
illustrated by the 1.1 reduced isotherm and the re- 
gions previously cited from the near-critical to the 
cirronoentherm. These difficulties will require a 
reassessment of the governing equations. 

As a further comparison, the solution of the 
governing equations for a second similar mixture 
(both ING), Mix? was carried out. The results are 
given in normalized coordinates in figure 6 . Com- 
paring the results for Mixi and Mix? indicates two 
possible similarities; ( 1 ) the slopes of the iso- 
therms and ( 2 ) the saturation locii are in better 
agreement than with the pure-pure background 
curves. Figure 9 gives a more direct comparison 
between Mixi and Mix? and corresponcing states. The 
6 * values for Mixi and Mix? were calculated using 
reference 13. At this time it does not appear rea- 
sonable to expect these curves to "collapse" to rep- 
resentative isotherms as composition dependence sig- 
nificantly alters the behavior near the saturation 
locus. Experimental two-phase choked flow data for 
fluid mixtures are indispensable. 

SUMMARY 

For the pure-pure problem involving fluids of 
simple (non-quantum) molecular structure, the con- 
ventional corresponding states approach to computa- 
tions of two-phase choked flows appears valid. The 
use of the van der Waals one fluid model ano a com- 
puter program such as Mollerup (ref. 10), with the 
appropriate reference fluio also appears to be a 
valid approach. 

The departures of two-phase choked flows in 
normalized coordinates for the pure-pure and the 
pure-mix problems are significant and composition 
dependent. 

For the mixture, the critical parameters, cir- 
condenbar and circondentherm neeo to be determined 
from the governing equations or experimental data. 
Once found, some normalization may be effective in 
grouping data as it becomes available. 
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TABLE 1. - MIXTURE COMPOSITION, NORMALIZING AND CRITICAL PARAMETERS FOR MIXj 


Composition-Mole Fraction 


1. Methane 0.9718 

2. Ethane 0.0477 

3. Propane 0.017 

4. 1 -Butane 0.0028 


5. n-Butane 0.0048 

6. 1 -Pentane 0.0013 

7. n-Pentane 0.001 

8. Hexane 0.001 


9. Heptane 0.0036 

10. Nitrogen 0.026 

11. Carbon dioxide 0,01 


Method 

V 

MPa 

K 

V 

cc/g 


m, 

g/g-mole 

G*. 

9 

g/cm -sec 

Circondenbar, 

MPa 

Circondentherm, 

K 

Eq. (16) 

4.6 

202.7 

5.79 

0.285 

18.02 

5280 

.... 


Ref. 13 

6.21 

213.0 

5.3 

0.335 

18.02 

5915 

— 


Mollerup code 
(ref. 10) 

6.46 

212.0 

4.26 

0.281 

18.02 

7340 

8.45 

267.4 


TABLE II. - MIXTURE COMPOSITION, NORMALIZING AND CRITICAL PARAMETERS FOR MIX 2 


Composition-Mole Fraction 

1. Methane 0.6975 3. Propane 0.092 5. n-Butane 0.029 

2. Ethane 0.156 4. 1-Butane 0.014 6. Nitrogen 0.0115 


Method 

p 

c* 

MPa 

V 

K 

V 

cc/g 

Z c 

m, 

g/g-mole 

g/cm-sec 

Circondenbar, 

MPa 

Circondentherm, 

K 

Eq. (16) 

4.53 

234.0 

5.36 

0.284 

22.75 

5455 

... 


Ref. 13 

8.95 

261.0 

4.39 

0.411 

22.75 

7040 



Mollerup code 
(ref. 10) 

8.9 

243.0 

3.09 

0.309 

22.75 

9650 

9.2 

299.0 


G* ■ /P c p c /Z c and GJ « |^h Gg for corresponding states applications with reference 
fluid (0) and working fluid (a) nearly equivalent for these fluids. 
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Figure 3. - Choked flow rate and pres- 
sure ratio ‘is a function of reduced 
pressure for a reduced temperature 
of 1.025 - conical nozzle. 
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Figure 5. * Primary mole fractions - pressure loci! at saturation conditions for mixture 
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Figure 6. - Temperature - entropy chart for mixture Mix*. See Table I for composition. 
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REDUCED INLET STAGNATION PRESSURE, P R#0 

Figure 7. - Reduced critical mass flux as a function of reduced inlet stagnation temperature 
for selected isotherms for mixture Mixj with conventional corresponding states locii as 
background. 
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REDUCED INLET STAGNATION PRESSURE. P R 0 

Figure 8. - Reduced critical mass flux as a function of reduced 
inlet stagnation pressure for selected isotherms for mixture Mix 2 . 

* B \ 
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1 

REDUCED PRESSURE, P R 0 

Figure 9. - A comparison of reduced critical mass flux as 
a function of reduced inlet stagnation pressure for se- 
lected isotherms for Mi*}, Mixo, and corresponding 
states using G* Tables I, II, ref. 13. 
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